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Abstract : The paper presents some thoughts about the plane strain problem of the 
viscous orthotropic composite materials cylinder under internal and external 
pressure with 

respect to using the average approximating method . To compute the interior stress , 
from the elastic solution we use the Volterra's principle and Sokolov's method in the 

corresponding integral equation to find the viscous solution . 



I. The axial symmetric plane strain problem of cylinder : 

We examine an orthotropic viscoelastic composite material cylinder which has the 
horizontal section within limit of 2 circles :r = a,r=b(a<b). 



z + 



x3 
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Choosing the cylindrical coordinates r , ^ , z ( the axial z is along with the 
cylinder) . The components of stress and deformation e -~ = E * '■ a * = °* 

are functions of r , t respectively . The two components of deformation-tensor : 



du(r 9 t) u{r,t) 

e r (r,t) = — ; £ (r,t) = 

or r 



(1.1) 



and the differential equation of equilibrium is : 

r da r (r,t) 

(1.2) 



dr 



+ <r r (r,t)-<T g (r,t) = 



when t = , boundary conditions : 
(1.3) 



. <r r (a,0) = -P ; a r (b,0) = -Q 



II . The Volterra Integral equation of the second kind : 

The displacement -differential equation of the cylinder in the case of viscoelastic 
plane-deformation : 



f & 



dr 



u(r, t) 



+ r\ Fr n(r,t) - 



Eu(r,t) 



(2.1) 



r -p \ 



E 

v < > J 



The general solution is : 
(2.2) 



u(r, t)'=C x r K +C 2 r 



•=0 



f T? \ 



E 

v v ' j 



c c 

1 , 2 are the arbitrary constants 



The elastic constants in (2.1) will be substituted by operators respectively , from the 
Volterra's principle we have : ([1]) 



2 d 2 u(r,t) du(r,t) E 



dr 2 



■+r- 



dr E 



-f-u(r 9 t) = 



(2.3) 

The equation (2.3) is rewritten as 
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-^- + a x (r) — -^- + a 2 (r)u(r,t) = 

dr dr 

/ x 1 ( \ E G m{t) 

a 1 (r 9 t) = - , a 2 {r,t) = — — — = — 

r ~ r £ r r 

(2.4) 



Assume that M ( a ^)- c o > n'Cfl.O-Q (j n [2] ) we obtain the solution of (2.4) 
from : 

(f>{x,t) + JK[(x,y);t]<f>(y,t)dy = f{x,t) 
(2.5) 



( x — v^ 
The kernel expression of : K[(x,y);t] = a l (x,t) + a 2 (x,t)- 



1! 



(2.6) 



And the formula of f(x,t) is : f(x,t)--C a l (x,t)-[c l x + C \i 2 (x,t) 

(2.7) 



1 m(x-y) 

K[(x, y);t] = 2 

In this case : * x 

(2.8) 



, lm (1/10). (2/5) 

( 100. + ? )t 

m := — 



100. +^/7 



(2.9) 



From the experimental test , if we have ° e ~ ° ' ^ (r ' 0) " Con ^ 

(2.10) 



Then : cr r (t) = e r .E r (t) 

(2.11) 
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From the results of (2.5) we obtain the analytical solution of u ( x > f ) j a e\ a ^) 
and ^(b.t) 

III. The approximate solution of the Volterra Integral equation of the second 
kind : 



3.1 The Sokolov's average approximate method 



The convergence speed of calculation process can be increased by Sokolov's 
method . The basic contents of this method is described as following : 



u(t) = v(t) + X \k{$, r)u(r)dr 
From the integral equation : 
(3.1) 



The n-order expression of u(t) is: u n (t) = v(t) + X JK(t,T){u n (r) + a n }dT 

(3.2) 

This relation has connection with the adjustment quantity : 

1 b 
a n =- \A n (t)dt 

b-a J a 

(3.3) 

The difference between two results «<■*)- «„ v) - u n _ x it) 

(3.4) 



a„ = 



2 b b 
■" ^O) \\ K{ t^n-^)-a n -Mrdt 
The recurrence relation ^ ' a a 

(3.5) 

Note that the convergence -condition of Sokolov's method is 



M = 



A(i-AKny l K(i-u) 



<i 



Here n is the project operator from the Banach 's space B into its subspace 
Bo ( u e B ) ([2]) 
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3.2 The flow chart of Sokolov's method 




Begin . 
Input the 
parameters. 




The initial solution 



Construction the integral 
equation 

b 

u(t) = v(t) + X ^K(t, z)u(z)dz 



T 



Find #1, define 

T 1 (x) = u 1 (x) 

Compute D = 



b b 

R(Z) = b-a-ZJJK(t 9 r)didt 



I 



define "«« , W 




i < N 



i = i + 1 




N 



Compute ^- l(x) ' ^ (x) 

and aM ' a ' +2 
estimate error sais °\. l \ 
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3.3. Problem definition - parameters and the analytical solution : 

In the following , modules ^(^0 and e\ r J) are given of the exponential 
expressions : 

a= 1 , b=2, C o=° • C 1= l. E e =0.5 , r = l . 

1. Parameters Information : 

:>m:=(l/x A 2)*(100.+(t/To) A (l/10))*(t/To) A (2/5)/(100.+(t/To) A (l/2));To:=l; 

We should choose n : = 4 " : = 5 " : = 6 



2. Activate the procedure : > xapxi ( m/x, -1 ,1 , 2 , 4 , K ) ; 



_ m(x-y) | 1 



* . To := 1 

5 

KETQUA- 



E a (r,t):= 



c_(x, := 



HOO A 
+ 1 



.0.1 



Ee 



ln(x)x(100. + / 1/10) )/ 2/5) l.(100. + / 1/10) )/ 2/5) x 



100. + Jt 



100. + ft 



0.002500000000 



+ 



0.1013313763 10 30 ( 100. + / 1/10) ) 12 / 24/5) 



( 100. + ft) 



12 



f 100 , ^ 



Eelx 



Comment : In this procedure note that we should choose the number n of 
recurrence from 4 up to 6 , so we get the result which is coincide to the solution 
of Schapery in 

[4] . In the comparison with the result in [3] the analytical solution of o e (a,i) and 
cr (b 9 t) obtained by Sokolov's method has a better smoothness . Figure 1 and 2 
describe 
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the convergence of s e (r 9 t) and e r (r 9 t) ■ Especially , if n>6 and using Maple 

version 9.5 it is a waste of time to treat the problem in detail , moreover the 
accuracy of 

solution will be influenced by accumulation of error in programme . 

Using Maple v ersion 6.0 the convergence of method can be found correctly when 
n = 5 . 
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8 




convergence of strain s r ( r >0 



Fig2 



Figure 3. and 4. show the manners of the graphs of stress coincide to the work 
of R.A.Schapery . To obtain these result , the most important work stage is the 

transformation differential equation to the correspondent integral equation and setting 
up the suitable boundary conditions. 
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graph of stress crJa,t) , a=1 



Fig 3 
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Fig4 



graph of stress a e( b ^) , £>= 7 



In Figure 5 we notice the schematic representation of the displacement at x = 1 

and x = 2 ; the manners of the graphs of strain s r ( r >0 and convergence 
process of 

s r ( r >0 according to the order { n -2 , n-1 ,n } can be shown in figure 6 . 
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10 
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F/g 5. graph of displacement u (r,t) 



Fig 6 graph of s r ( r »0 . 
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Fig 7 . graph of £ & ( r ^ 



Fig 8 . graph of ^(r.O 
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12 



IV . Conclusion : 

The plane strain problem of the viscous orthotropic composite materials cylinder 
under internal and external pressure can be solved by other approximate methods : 
direct , collocation , quasi-elastic , Laplace transform inversion 

. . . ([3]) . Almost the algorithm regularly depends on the displacement differential 
equation and the establishment of procedure are also carried out from it . The 
Average Approximating Method on Functional Adjustment 

Quantity ( Sokolov's method ) and the Volterra's principle corresponding to the 
integral equation may be used to find the viscous solution . Moreover this 
application makes increasing for the convergence speed of the 

solution u n it) . From the first approximation of the solution u , we find the 
adjustment quantity for the next and so on . Note that we should choose the number 
n of recurrence and the expressions of a 1 (x,t) , a 2 (x 9 t) 

appropriately to get the analytical result which is satisfied the given boundary 
conditions of the problem . 
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